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Abstract 


In  this  paper,  we  are  concerned  with  the  per  packet  transmission  delays 
induced  by  the  Capetanakis  collision  resolution  protocol  for  infinitely  large 
number  of  identical  bursty  users. 

We  first  correct  the  existing  lower  and  upper  bounds  on  the  expected  per 
packet  transmission  delay  E{d}. 

Then,  we  proceed  by  developing  a  new  upper  bound  on  E{d}.  This  new  bound 
is  much  tighter,  approaching  the  lower  bound  for  large  arrival  rates. 
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Introduction 


In  this  paper,  we  are  concerned  with  the  "multiple-access"  problem, where  a 
large  number  of  independent,  packet  transmitting,  bursty  users  request  access  to 
a  common  channel.  We  consider  "random-access"  transmission  protocols,  as  more 
efficient  for  the  present  problem  [4]. 

In  more  specific  terms,  the  user  and  channel  models  considered,  are  as 
follows: 

(i)  The  users  are  independent  from  each  other,  and  they  can  communicate 
with  each  other  only  through  the  channel.  Furthermore,  the  users 

are  large  in  number,  they  are  identical  and  bursty,  and  each  generates 
packets  of  fixed  common  length.  The  cumulative  input  to  the  channel 
traffic  is  a  Poisson  process. 

(ii)  The  common  channel  is  perfect,  i.e.  there  are  no  channel  errors.  In 
addition,  the  channel  time  is  divided  in  slots  of  Identical  length, 
where  this  length  is  equal  to  the  length  of  one  packet. 

(iii)  The  feedback  channel  is  perfect,  i.e.  it  does  not  induce  propagation 
delays.  Furthermore,  the  feedback  channel  broadcasts  with  no  errors 
the  outcome  from  each  channel  slot.  In  particular,  it  broadcasts  a 
trenary  sequence,  where  the  value  of  each  digit  from  this  sequence 
indicates  if  the  corresponding  slot  was  empty,  busy  with  exactly  one 
packet,  or  busy  with  at  least  two  packets. 

(iv)  The  transmission  characteristics  imply  low-level  synchronization  among 
the  users.  In  particular,  each  user  is  allowed  to  transmit  at  most 
one  packet  at  the  time,  and  he  can  start  transmission  only  at  the 
beginning  of  some  channel  slot.  Thus,  some  channel  slot  is  empty  if 
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no  user  transmitted  a  packet  within  it.  Some  channel  slot  is  busy 
with  exactly  one  packet  if  exactly  one  user  transmitted  within  it. 

*  It  is  assumed  that  xn  this  last  case  the  transmitted  packet  is  received 
correctly.  Incorrect  transmission  occurs  only  if  at  least  two  packets 
are  transmitted  within  the  same  slot.  Then,  collision  occurs,  and  the 
involved  packets  are  lost  completely.  If  so,  those  packets  are  restored 
in  the  queues  of  the  corresponding  users,  and  are  retransmitted  within 
some  future  slots. 

Any  "random-access"  transmission  protocols,  which  are  appropriate  for  the  above 
model,  are  characterized  by  a  number  of  performance  parameters.  The  number  one  such 
performance  parameter,  which  in  fact  determines  the  eligibility  of  a  given  random- 
access  transmission  protocol,  is  stability.  A  stable  random-access  transmission 
protocol  maintains  the  rate  of  the  cumulative  input  Poisson  traffic,  while  an 
unstable  such  protocol  does  not.  Given  a  stable  random-access  transmission  protocol, 
two  other  performance  characterisitcs  for  comparison  with  other  stable  such  protocols, 

are  the  throughput  and  the  per  packet  transmission  delay.  The  throughput  is  defined 

..  //successful  transmissions  ......  ,  _  . 

as  the  ratio  p  =  -vrr - : - : - — - r -  induced  by  the  random-access  transmission 

//transmission  attempts 

protocol,  and  it  is  related  to  the  channel  capacity.  In  fact,  the  channel  capacity 
is  the  maximum  of  all  throughputs  induced  by  stable  random-access  transmission 
protocols.  The  per  packet  transmission  delay  is  defined  as  the  time  between  the 
arrival  of  some  packet  and  its  successful  transmission,  where  time  is  measured  in 
number  of  channel  slots. 

The  oldest  existing  random-access  transmission  protocol  for  the  user  and  channel 
models  stated  in  this  paper,  is  the  slotted-Aloha.  The  problems  regarding  the  opera¬ 
tion  of  the  slotted-Aloha,  are  by  now  well  known  [1-4],  and  we  will  not  discuss  them 
here.  Instead,  we  will  focus  on  the  protocol  by  Capetanakis  [1-4].  Capetanakis' 
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protocol  has  been  called  collision-resolution  protocol,  it  is  stable  for  input 
Poisson  rates  below  .3465,  and  it  induces  a  throughput  of  .43.  Refinements  of 
Capetanakis'  algorithm  by  Gallager,  Massey,  and  Mosely,  have  increased  the  input 
Poisson  rate  for  stability  to  .375,  and  the  induced  throughput  to  .488  [4]. 

In  this  paper,  we  concentrate  on  the  per  packet  transmission  delay  induced  by 
Capetanakis'  collision  resolution  protocol.  Afterall,  the  very  reason  for  considera¬ 
tion  of  random-access  protocols  is  the  improvement  of  delays. 

Capetanakis  studied  the  per  packet  transmission  delay  (or  waiting  time)  induced 
by  his  otocdL.  The  resulting  expression  being  complex,  a  lower  and  an  upper  bound 
on  this  delay  were  developed.  The  procedure  for  the  derivation  of  these  bounds  can 
be  found  in  [1,4]. 

Here,  we  develop  tighter  lower  and  upper  bounds  on  the  per  packet  transmission 
delays  induced  by  Capetanakis’  protocol,  through  the  correction  of  some  step  in  the 
used  procedure  [1,4]  first,  and  then  through  a  fresh  approach. 


The  Capetanakis  Protocol  -  Notation 

We  will  use  reference  [4],  since  we  feel  that  Capetanakis'  protocol  and  its 


analysis  are  best  explained  there.  We  will  use  basically  the  same  notation  as  in 


[4],  and  we  will  describe  the  Capetanakis  Collision  Resolution  Protocol  (CCRA)  only 
briefly  and  quantitavely. 


The  CCRA  is  activated  just  after  a  collision  slot.  Then,  through  the  feedback 
channel,  all  users  are  instructed  to  withhold  newly  generated  packets  until  the 
collision  is  resolved.  The  collision  resolution  interval  (CRI)  is  the  number  of 


slots  required  for  the  resolution  of  the  collision,  that  is  the  successful  trans¬ 


mission  of  all  packets  involved  in  the  collision.  The  collision  resolution  is 


obtained  through  the  application  of  the  CCRA,  during  the  CRI.  The  CCRA  is  based  on 
the  following  general  principle: 
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After  each  collision  slot  within  the  CRI,  each  user  involved 
in  this  collision  flips  independently  a  binary  fair  coin  with 
outcomes  0  and  1.  Among  those  users,  only  the  ones  with  out¬ 
come  0  transmit  within  the  next  slot.  Until  the  initial 
collision  is  completely  resolved,  the  resolution  of  no  other 
collision  within  the  CRI  is  attempted. 

The  CCRA  has  the  structure  of  a  binary  tree,  where  each  node  within  the  tree 
is  taken  to  the  leaves-depth,  before  other  nodes  on  the  same  depth  with  the  original 
node  are  resolved. 

Using  the  reasonable  for  infinite  number  of  identical  bursty  users  assumption, 
that  during  any  CRI  there  may  be  at  most  one  new  packet  arrival  per  user,  we  define 
as  in  [4]  the  following  parameters: 

:  The  expected  length  of  a  CRI,  given  that  the  number  of  packets  involved 
in  the  collision  within  the  slot  with  which  the  CRI  starts,  is  equal  to  N. 
:  The  number  of  packets  involved  in  the  collision  of  the  first  slot  of  the 
ith  CRI,  from  the  beginning  in  time  that  the  system  starts  operating. 

:  The  length  of  the  ith  CRI. 

X^  :  The  X^  for  i-*00,  that  is  in  steady-state. 

Y 

00  :  The  Y.  for  i-*». 
x 

Y  :  In  steady-state,  the  length  of  the  CRI  in  progress,  when  some  new  packet 
fl 

arrival  occurs. 

Yj  :  In  steady-state,  the  length  of  the  CRI  in  progress,  when  some  packet 
departs  from  the  system  (that  is,  when  the  packet  is  successfully 
transmitted) . 

X,  :  In  steady-state,  the  number  of  packets  involved  in  the  collision  of  the 

d 

first  slot  of  the  CRI,  during  which  a  packet  departs  from  the  system. 
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D  :  In  steady-state,  the  time  interval  measured  in  number  of  slots,  from  the 
arrival  of  some  packet  to  its  departure  from  the  system. 

X  :  The  rate  of  the  cumulative  input  Poisson  process. 

1  ;  if  i  =  N 

V5 

0  ;  otherwise 

We  will  point  out  here  that  the  existence  of  steady-state  for  the  CCRA  has  been 
proven  rigorously  by  Capetanakis. 

Using  the  notation  presented  in  the  present  section,  in  the  next  section  we 
will  outline  the  approach  taken  and  the  results  obtained  in  [4] ,  regarding  the  per 
packet  transmission  delays  induced  by  the  CCRA. 

3.  The  CCRA  Transmission  Delays 

In  [4],  lower  and  upper  bounds  on  the  expected  per  packet  transmission  delays 
(waiting  times)  E{d}  induced  by  the  CCRA,  have  been  obtained.  The  bounds  are  valid 
for  rates  X  in  the  region  [0,  .3465].  Within  this  X  region,  the  expected  length 
of  each  CRI  is  finite,  and  the  system  reaches  steady-state. 

The  following  approach  was  taken  in  the  development  of  the  bounds  (in  [4]). 

First  a  recursive  expression  for  the  evaluation  of  was  obtained.  Based  on 
this  expression,  the  following  lower  and  upper  bounds  on  were  obtained: 

Ln  >  2.8810  N  -  1  +  2  6qn  -  0.8810  (1) 

LN<  2.8867  N+  6qn  -  1.8867  6lN  (2) 

In  parallel,  it  was  also  found  that  the  following  equality  holds: 

E{Y^} 

E{V  -  eTTJ 


(3) 
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From  this  point  on,  the  key  expression  used  for  the  eventual  evaluation  of 
the  bounds  for  E{d},  was  the  following: 


I 'A -4 


Expression  (4)  is  true  due  to  the  fact  that  the  arrival  process  is  Poisson, 
and  it  results  in  the  following  expression: 


E{X , }  =  X  E{Y  } 
d  a 


The  bounds  in  (1)  and  (2),  in  conduction  with  expression  (5)  and  the  obvious 

fY,  /  ) 


relation  E 


X,  =  N 
d 


=  L^,  result  then  in  the  following  lower  and  upper  bounds  for 


the  expected  value  ECy^}: 


E{Y, }  >  2.8810  X  E{Y  }  -  1.8810  +  2.8810  e 
d  —  a 


-X  E{Y  } 


E{Y,}  <  2.8867  X  E{Y  }  +  (1  -  1.8867X)  e“  '  (7) 

d  —  a 

The  bounds  in  (6)  and  (7)  are  expressed  as  functions  of  the  rate  X  and  the 
expected  value  E{Y }. 

Subsequently  the  following  lower  and  upper  bounds  on  the  expected  delay  E{d} 

were  found,  in  terms  of  E{Y  }  and  E(Y,}: 

a  d 

4  E{Y  }  +  4  E{Y,  -  1}  <  E{d}  <  4  e(Y  }  +  E(y,  -  1}  (8) 

2a2d  —  —  2a  a 

The  bounds  in  (8) ,  in  conjuction  with  the  bounds  in  (6)  and  (7) ,  result  clearly 


in  the  following  bounds  for  e{d}  : 


e(d)  >  (|  +  1.4405  X)  E{Y, 


/  -X  e{y  } \ 
}  -  1.4405  \  1  -  e  / 


e(d}  <  (4  +  2.8867  X)  E{Y  }  +  (1  -  1.8867X)  e~X  -  1 

—  2  a 

t 

In  expressions  (4.42)  and  (4.43)  in  [4],  a  parenthesis  is  missing. 
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The  bounds  in  (9)  and  (10)  are  functions  of  A  and  E{Y  }  only. 

3 

2 

Finally,  tight  bounds  on  EfY^}  and  e{Yoo),  result  in  tight  lower  and  upper 
bounds  for  e{y  }  in  (3) .  These  bounds  clearly  provide  lower  and  upper  bounds  for 

3 

E{d}  through  (9)  and  (10),  and  they  are  given  by  the  following  expressions: 


E{Y  }  < 
a  — 


5.964A  -f  1  -  5.964A(1-A) (1  -  2.8867A)/fl  -  2.8867A(1-A) 1 

1  -  8.333A2 


5. 964 A  +  1 
1  -  8.333A2 


.22  <  A  <  .3464 


A  <  .22 

(ID 


I  5.897A  +  (1  -  5.897A)  (1  -  2.8867A)/fl  -  2.8867A(1-A)  1  . 

- - -  ;  A  <_  .1696 

1  -  8.300AZ 

(12) 

- - - r-  ;  .1696  <  A  <  .3464 

1  -  8 . 300A^ 


The  bounds  in  (11)  and  (12) ,  in  conjuction  with  the  bounds  in  (9)  and  (10)  , 
provide  lower  and  upper  bounds  for  E{d},  which  are  functions  of  the  rate  A  only, 
as  desired. 

These  bounds  have  been  computed  for  different  A  values  [4].  We  will  discuss 
the  results  from  the  computations  in  [4]  as  they  compare  to  our  results,  in  the 
following  section. 

4.  Modified  E{p}  Bounds 

In  section  3,  we  outlined  the  approach  taken  in  [4]  for  the  computation  of 
lower  and  upper  bounds  on  the  expected  per  packet  transmission  delay  E{d}. 

The  key  expression  for  the  computation  of  the  E(d}  bounds  in  [4],  was  expression 
(4)  in  section  3  of  the  present  paper.  But  the  assumption  throughout  the  whole 
related  analysis,  is  that  there  exists  at  least  one  packet  arrival  within  the  Ya 
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in  (4). 


slot  period.  Thus,  the  relevant  key  expression  in  the  computation  of  the  E{d} 

bounds  should  be  an  expression  for  the  expectation  E{x, |y  =  L,  X,  >  l},  rather 

(Xd/  }  d|  *  d" 

than  the  expression  for  the  expectation  E<  L  _  in  (4). 

In  this  section,  we  will  substitute  expression  (4)  by  an  expression  for  the 

jxd  /  { 

expectation  Es  /  ,  an^  ue  subsequently  compute  modified  lower 

(  /  Ya  =  L’  Xd  -  M 

and  upper  bounds  for  the  expected  value  E{d} . 

Due  to  the  fact  that  the  packet  arrival  process  is  Poisson  with  rate  X,  we 


clearly  have: 


X ,  =  k 
a 


Y  =  L,  X,  >  1 
a  d  — 


0  ;  k  =  0 


J  ; 

(■■ ;  y, .  j  '  “ 


-XL  (XL) 


;  k  >  0 


From  expression  (13),  we  have  in  a  straight-forward  manner: 


Ya-1,  X(J>1)  1-e 


j 


-XY 
e  a 


(XY  V 
a 


.  -XY 
1-e  a 


X ,  >  1 
d  — 


;  where  the  expectation  in  (15)  is  with  respect  to  Y  . 

3L 

At  this  point,  let  us  observe  that  the  implicit  assumption  is  that  at  least 
one  arrival  occurred  during  the  slot  interval  represented  by  Y  .  This  assumption 
implies  X^  1.  It  is  due  to  this  observation  that  Massey  [4]  dropped  the  condi- 
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tioning  X^  ^  1  from  the  expectation  E 


IV,  -  ^  ■!• 


Unfortunately,  however. 


the  above  observation  does  not  imply  that  the  X^  process  is  Poisson.  The  distribu¬ 
tion  of  X^  is  rather  reflected  by  expression  (14) . 

Using  the  above  observation,  we  will  delete  from  now  on  the  conditioning 
X^  >  1  from  expressions  involving  expectations  on  Y^.  Then,  the  X^  1  condi¬ 
tioning  in  (15)  is  deleted,  and  we  proceed  as  follows. 

From  expression  (15)  we  obtain  directly  the  following  expressions: 


(16) 


(17) 


In  expressions  (16)  and  (17)  the  expectations  are  with  respect  to  Y  . 

3 

Now,  we  will  use  the  upper  bound  for  L  ,  as  given  by  expression  (2),  where  the 

ft 

term  6^  is  deleted  due  to  the  X^  ^  1  conditioning.  We  obtain  then,  the  following 
inequality 

|  =  Ln  <  2.8867  N  -  1.8867  5lfJ  ;  N  >  1  (18) 

Directly  from  (18),  we  also  obtain  the  following  expression. 

<i9> 


Substituting  expressions  (16)  and  (17)  in  expression  (19) ,  we  obtain: 


AY 


Xd>  u<  2.8867  E( 


1  -  e"XYa 


-  1.8867  E 


iw  I 


(20) 
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In  the  right  part  of  inequality  (20) ,  we  add  and  subtract  the  term  E 
We  then  obtain: 


v  .  .  ,  <  2.8867  X  E(Y  }  +  E 
xd  >  1 1  -  a 


XY 


XY 

e  a-l 


w  -XY 
XY  e  a 
a _ 

i  _Xy 

1  -  e  a 


(21) 


Observing  now  that  Y  >  1,  and  that  for  x  >  0  the  function 


Xx 


Xx  . 

e  -1 


is  monotoni- 


cally  decreasing  with  increasing  x,  and  it  is  convex  for  x  _>  1,  we  can  obtain  the 
following  bound: 


XY. 


XY 


\< 


A  7 

e  -1 


;  For  Y  >1 
a  — 


(22) 


Substituting  the  bound  in  (22)  in  expression  (21),  we  obtain: 


Expression  (23)  provides  an  upper  bound  on  the  expectation  E 


1  ‘d 


I 


(23) 


which  is  different  than  the  upper  bound  in  (7),  if  the  conditioning  X,  >  1  is 

d  — 

deleted  as  necessary  and  implicit. 

Using  the  lower  bound  for  L  ,  as  given  by  expression  (1),  we  can  derive  similarly 

Yd  /  ) 
a  lower  bound  fdr  E  {  .  ' 

X ,  >  l 
a 


I  /xdH) 


Indeed,  deleting  again  the  term  in  (1),  we  obtain: 


I 


y 


Xd-N,  xd  >  1  -  h  i  2-8810  »  ‘  1  -  °-8810  S1N 


(24) 


Averaging  out  with  respect  to  N,  as  in  (19),  and  substituting  expressions  (16) 
and  (17),  we  obtain  from  (24): 


I 


>  2.8810  E 


I.  XY 


!l  -  <rXYa  \ 


h- 


0.8810  E 


,v  -XY 
AY  e  a 
a _ 


1  -  e  a 


(25) 
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•  w  -AY  . 

XY  e  as 

Adding  and  subtracting  the  term  2E  { — - - >  in  the  right  part  of  (25),  we 

-XY.  J 


1  -  e  a 


obtain: 


Jv 

I 


>  2.8810  A  E{Y  }  -  1  +  2  E 


iY.  I 


TT 


71 


(26) 


Due  to  the  convexity  of  the  function 


Ax 


Ax  , 
e  -1 


for  x  >  1,  and  due  to  the  fact  that 


Y  >  1,  applying  the  Jensen  inequality  we  obtain  the  following  bound: 


.K  I  > 
1”*-!  f  - 


e{y  } 

i  a 

7e(y  } 

a 


(27) 


Substituting  the  bound  in  (27)  to  the  inequality  in  (26),  we  finally  obtain: 


(28) 


(  Y  /  ) 

Expression  (28)  provides  a  lower  bound  on  the  expectation  * 

which  is  different  than  the  lower  bound  in  (6),  if  the  conditioning  >  1  is 
deleted  as  necessary  and  implicit. 

Using  now  the  bounds  in  (23)  and  (28),  in  expression  (8),  after  deleting  the 
conditioning  >_  1,  we  obtain  the  following  modified  bounds  for  the  expectation 
E{d)  : 


,  Ae(Y  } 

E(D}  >  (7  +  1.4405  A)  E{Ya)-l+-Wf- 


-1 


E{D}  <  4  +  2.8867  A)  e{y  }  -  1  +  — r- 
—  t  a  A 


e  -1 


(29) 


(30) 


The  bounds  in  (29)  and  (30)  are  different  than  the  corresponding  bounds  in 
(9)  and  (10),  and  they  are  functions  of  the  rate  A  and  the  expected  value  e(Y  }. 
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If  the  bounds  on  the  expected  value  E{Y  },  as  given  by  expressions  (11)  and  (12), 
are  used  in  expressions  (29)  and  (30),  modified  lower  and  upper  bounds  on  E{d} 
can  be  obtained.  These  last  bounds  are  functions  of  the  rate  X  only,  as  desired. 

We  performed  parallel  computations  for  the  bounds  in  (9)  and  (10)  and  the  bounds 
in  (29)  and  (30),  for  different  values  of  X.  Our  results  are  exhibited  in  table  1. 

As  expected,  due  to  the  correction  we  obtained  by  the  addition  of  the  conditioning 
Xj  1  1.  °i*r  modified  lower  and  upper  bounds  are  higher  than  the  corresponding  bounds 
in  [4],  for  all  X  values.  Also,  as  expected,  the  correction  has  stronger  effects  on 
low  X  values.  In  fact,  for  X  values  below  0.15,  our  modified  lower  bound  is  higher 
than  the  upper  bound  in  [4].  In  general,  the  curves  indicating  the  modified  lower 
and  upper  bounds  are  shifted  upward  versions  of  the  corresponding  curves  in  figure 
4.1  of  [4]. 


X 

Lower  Bounds 

Upper  Bounds 

0 

M 

0 

M 

0.0500 

0.521 

0.569 

0.531 

0.645 

0.1000 

0.599 

0.701 

0.664 

0.881 

0.1500 

0.780 

0.961 

1.009 

1.319 

0.1696 

0.902 

1.125 

1.252 

1.595 

0.2000 

1.200 

1.505 

1.842 

2.236 

0.2500 

3.423 

3.958 

5.762 

6.231 

0.3000 

9.009 

9.521 

14.560 

15.096 

0.3333 

38.067 

38.508 

58.204 

58.781 

0.3400 

80.203 

80.643 

121.458 

122.043 

0.3450 

371.910 

372.351 

559.442 

560.032 

0  :  Bounds  in  [4] 

M  :  Modified  Bounds  in  (29)  and  (30). 


Table  1 
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5.  Improved  Bounds  for  E{d) 

In  section  4,  we  simply  corrected  the  bounds  In  [4]  for  E{d}. 

In  this  section,  we  will  take  a  fresh  approach,  to  develop  tighter  bounds  on 


We  will  first  develop  an  exact  expression  for  a  per  packet  waiting  time  para¬ 
meter.  We  will  then  present  an  intuitive  analysis  and  subsequent  bounds  on  this 
parameter.  Finally,  using  these  bounds,  we  will  compute  bounds  on  E{d}. 

Given  that  N  packets  are  involved  in  the  collision  within  the  first  slot  of 
some  CRI,  we  will  first  seek  expressions  for  the  expected  transmission  delay  for 
each  of  the  N  packets,  within  the  CRI.  We  will  adopt  Capetanakis'  collision  resolu¬ 
tion  protocol  [1  -  4]. 

After  the  first  collision  slot  within  the  CRI,  which  started  with  N  collided 
packets,  let  i  out  of  the  N  users  have  outcome  0  from  flipping  their  fair  coin. 


Then,  the  corresponding  outcome  for  the  remaining  N-i  users  is  1.  Let  us  denote 

this  event  (i,  N-i).  Then,  given  the  event  (i,  N-i),  each  of  the  N  packets  gives 

i  N-i 

outcome  0,  with  probability  — .  It  gives  outcome  1,  with  probability 

Define: 

d*  :  The  expected  delay  in  the  transmission  of  some  packet  within  the 
;j  =  0,1  CRI,  given  that  the  packet  gave  outcome  j  at  the  first  trial,  and 
that  the  first  trial  event  is  (i,  N-i). 


“N/i 


"N 


The  expected  delay  in  transmission  of  some  packet  within  the  CRI, 
given  that  N  packets  were  involved  in  the  collision  within  the  first 
slot  of  the  CRI,  and  that  the  first  trial  event  is  (i,  N-i). 

The  expected  delay  in  transmission  of  some  packet  within  the  CRI, 
given  that  N  packets  were  involved  in  the  collision  within  the  first 


slot  of  the  CRI. 
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Then,  we  have: 


.i  i  ,  ,1  N-i 

mN/i  d0  N  +  dl  N 


(31) 


Also,  defining  =  0,  observing  that  tn^  =  0,  and  considering  the  operation 
of  the  CCRA,  we  can  easily  see  that  the  following  expressions  hold  for  N  >  2: 


d0  "  1  +  “i 


dl  =  i  +  4  +  V. 


(32) 


;  where  is  defined  in  section  2  of  this  paper.  Substituting  expressions  (32)  in 
(31)  we  obtain: 

Vt  '  1  +  “1  5  +  T  Vi  +  V  •  "1  2  (33) 

Clearly,  the  probability  with  which  the  event  (i,  N-i)  occurs  is  (*?)  -i— . 

2N 

Therefore,  we  have: 


*N  (i)  b 


“N/i 


i=0 


Substituting  expression  (33)  in  expression  (34),  we  obtain: 


(34) 


N-l 


N-l 


(2 


i=l 


i=0 


(35) 


Now,  we  use  the  following  expressions  from  [4): 

N 


2"  L«  ■  2”+  2  (  i)  !  K-2 


i=0 


Lo  ■  H  ■  1 


(36) 
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Using  expressions  (36)  and  (35),  we  finally  obtain  the  following  exact  recursive 
expressions  for  m^,  which  also  involve  the  L^'s: 


(37) 


In  the  appendix,  we  present  an  alternative  procedure  for  the  exact  evaluation 
of  the  m^  parameters. 

Unfortunately,  the  expressions  in  (37)  are  complex,  and  the  development  of 
bounds  for  the  m^'s  from  them,  seems  impossible.  For  that  reason,  we  will  use  some 
intuitive  analysis  for  the  evaluation  of  the  m^'s.  Then,  we  will  verify  the  results 
from  the  intuitive  analysis,  by  computationally  comparing  with  the  exact  expressions 
in  (37). 

a.  Intuitive  Analysis  for 

Consider  some  CRI,  which  starts  with  a  slot  with  2N  collided  packets.  If 
N-*”,  then  with  probability  close  to  one,  the  first  trial  will  result  in  N  outcomes 
equal  to  0,  and  N  outcomes  equal  to  1.  Thus, 


m2N 

N-*®° 


1  +  m^  ;  with  probability  ~  1_ 

N  2 


(38) 


1  +  +  m^  ;  with  probability  - 

Therefore,  we  conclude  from  (38)  that  for  N-*00  we  have: 


m2N  ~  1  +  ®N  +  LN/2 


(39) 
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Following  the  recursions  induced  by  expression  (39),  we  obtain: 


k-1 


+  +  i  ^  ]  L  £  ;  for  N  power  of  2,  and  N  large 


(40) 


i=l 


Similarly,  for  N  equal  to  a  power  of  2,  and  N  large  enough,  we  also  obtain: 

k-1 


Li  +£  \i 


L  ,  -  k  +  L, 

k  1  ^ 

Z  i=l 


(41) 


Substituting  (41)  in  (40),  we  then  obtain: 


It  L1 

m2k  "  2  +  ~T  "  T  +  mi 


(42) 


if 

Denoting  N  =  2  and  k  =  £og2N,  we  finally  obtain  from  (42)  the  following 
approximation: 

*  .5  £og2N  +  ~  +  C  ;  N  -*>  » 


(43) 


;  where  C  is  a  variable  equal  to  -  y-  ,  whose  influence  on  the  expression  (43) 
is  infinitesimal  for  large  N's. 

Let  us  now  define: 


m^  =  .5  £og2N  +  y— 


(44) 


Then,  ra^  provides  an  approximation  for  the  parameter  m^,  whose  exact  expression 
is  given  by  (37) . 

It  is  unknown,  at  this  point,  how  good  an  approximation  of  n^,  m^  is.  We  will 
study  the  closeness  between  m^  and  ra^  numerically.  We  present  our  results  on  that, 
in  the  following  subsection. 
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V 

r 


a 


b.  Comparison  between  m^  and  m^ 

We  denote  by  m^  the  exact  expression  in  (37) .  We  denote  by  the  approxima¬ 
tion  in  (44).  We  define: 


'  "k  "  >  (45) 

We  study  the  closeness  between  the  expressions  for  m^  and  numerically,  by 
computing  the  parameters  e^  and  N  ^e^  for  different  N’s.  The  results  of  our  compu¬ 
tations  are  exhibited  in  table  2. 

From  the  results  in  table  2,  we  observe  that  for  increasing  N,  the  parameter 
eN  increases  monotonically,  while  the  rate  of  increase  N  ^eN  decreases  monotonically . 
In  fact,  we  observe  that  for  N  >_  2,  the  N  ^e^  values  remain  within  the  interval 
[0,  .019].  Due  to  the  monotonicity  of  the  function  f(N)  =  N  ^e^,  the  observed 
upper  bound  .019  on  f(N)  is  pessimistic  for  N  values  larger  than  70.  Using  this 
pessimistic  bound  for  large  N's, together  with  the  lower  bound  0  for  f(N)  and  with 
expressions  (44)  and  (45) ,  we  obtain  the  following  lower  and  upper  bounds  on  the 
exact  n^,  for  N  1: 

L  L 

.5  £og2N  +  j-  -  .5  <51n  <  <  .5  £og2N  +  ~  +  .019N  -  .5  «1N  (46) 

The  bounds  in  (46)  will  be  used  for  the  computation  of  lower  and  upper  bounds 
for  the  expectation  E{d}. 

c.  Lower  and  Upper  Bounds  on  E{d} 

We  will  first  proceed  with  bounds  for  the  expectation  Eim^}.  These  bounds  will 
be  found  through  the  computation  of  the  expected  values  of  the  bounds  in  (46) ,  condi¬ 
tioned  on  >  1  (as  in  section  4  of  this  paper) . 

We  will  first  compute  the  expected  values  of  each  term  in  the  bounds  appearing 
in  (46).  We  recall  that  the  ramdon  variable  N  corresponds  to  X^  =  N,  where  X^  is 
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N 

% 

a 

“n 

1 

0.0  .  " 

0. 

5C00GCC+00 

2 

0 .3000000+0  1 

0. 

30000CC+01 

3 

0.466667D+01 

0. 

46253 1C+01 

4 

0.6333330+01 

0. 

6261900+01 

5 

0. 7961900+01  -- 

!  0. 

,7870490+01 

6 

0.9’55361D+01 

i...  o. 

9449C1D+01 

'7 

‘  9. 1111790+02 

:  o. 

1 1004 1D+02 

S 

0.1266310+02 

0. 

1254270+02 

9 

0.1-4 19500  +  02 

0. 

1406950+02 

10 

0.1-571730+02 

0. 

1558760+02 

11 

0.1723190+02 

0. 

,  1709880  +  02 

12 

0. 1574020+02 

0. 

,1360440+02 

13 

0.2024330+02 

0. 

2C105CC+02 

14 

0.2174170+02 

0. 

2160140+02 

15 

0 . 2  3236  ID  +  02 

0. 

2309410+02 

16 

0.247?70D+02 

0. 

,2 453340+02 

17 

0.2521480+02 

0. 

26C693C+0 2 

13 

0.2769990+02 

0. 

,  2755370  +  02 

19 

0.2918250+02 

0. 

.2903530+02 

20 

0.3066310+02 

0. 

.3051490+02 

21 

0. 3214170+02 

0. 

3199270+02 

22 

0.3361370+02 

0. 

,3346890+02 

23 

0.3509410+02 

0. 

,3494260+02 

24 

0.3556820+02 

0. 

,3641700+02 

25 

0.3504100+02 

0. 

,3768910+02 

26 

0.3951250+02 

0. 

,3936010+02 

27 

0.4098300+02 

0, 

,4083010+02 

28 

0.4245250+02 

0. 

,4229910+02 

29 

0.4392100+02 

0. 

,4376720+02 

30 

0.4533870+02 

0. 

,4523440+02 

31 

0.4555540+02 

0, 

,4670060+02 

32 

0.4532140+02 

0, 

,  4816640+02 

33 

0.4973670+02 

0  i 

,4963120+02 

34 

0.5125130+02 

0. 

,5109550+02 

35 

0.5271520+02 

0. 

.5255910+02 

36 

0.5417650+02 

0, 

,5402210+02 

37 

0.5564120+02 

0, 

,5545460+02 

33 

0.5710330+02 

0, 

.5594650+02 

39 

Q. 5356500+02 

0, 

.5840790+02 

40 

0.6002610+02 

0, 

.5986580+02 

41 

0.6 146680+32 

0, 

.6132930+02 

42 

0.6294710+02 

0, 

.6278930+02 

43 

0.6440690+02 

0, 

.6424390+0 2 

44 

0.65S663D+02 

0, 

.6570320+02 

45 

0.6732540+02 

0, 

.6716710+02 

46 

0.6378410+02 

0 

.6852560+02 

47 

0. 7G2425D+02 

0 

.7008380+02 

48 

0.7170050+02 

0, 

.7154170+02 

49 

0.7315820+02 

0 

.7299920+02 

50 

0. 7461560  +  02 

0 

.7445650+02 

51 

0.7507280+02 

0 

.7591350+02 

52 

0.7752960+02 

0 

,7737020+02 

53 

0.739862D+02 

0 

.7832570+02 

54 

0.6044250+02 

0 

.8025290+02 

55 

0.8189360+02 

0 

.8173=80+02 

56 

0. 6335450+02 

0 

.6319460+02 

57 

0,8 48 1010+02 

0 

.8455000+02 

58 

0.8526540+02 

0 

.6610530+02 

59 

0. 6772060+02 

0 

.8756040+02 

60 

0.891 7560+02 

0 

.8901520+02 

61 

0.9  36303D  +  02  : 

0 

.  SC  46990  +  92 

62 

0.9203460+02 

0 

.9192430+02 

63 

0.9353920+02 

0 

.9337560+02 

64 

0.9499340+02 

0 

.9433260+02 

65 

0. 9544730+02 

0 

.9623550+02 

66 

0.97901 1.0  +  02 

0 

.9774030+02 

67 

0.9935460+02 

0 

.9919380+02 

68 

0. 1006080+03 

0 

.  100647C  +  03 

69 

0.1022620+03 

0 

.  102 1CCC  +  03 

70 

0.1037150+03 

0 

.  1035530  +  03 

MAXNO  =  5  MAX  =  .182834D  -  01 

MINO  *  2  MIN  =  .111022D  -  15 


-0.5000000+00' 
0.222045D-15  . 
0.136174D-01 
.  0.17857 1  D-Ol 
;  0.1 8  22  34D“0 1 
>  0.174334D-01 
O.162474D-0I 
0.15Q5310-01 
0.139529D-01 
0. 1296910-01 
0.120967D-01 
0.1132330-01 
0. 1063560-01 
0.1002170-01 
0.9471 110-02 
0.897530D-02 
0.8526910-02 
0.8119820-02 
0 .7748820-02 
0.7409510-02 
0.7098120-02 
0  .68 1 1 43Q-02 
0. 6546690-02 
0.630152  D~C2 
0.6073840-02 
0. 5861880-02 
0  .5664  08  0-02 
0.5479070-02 
0.5305670-02 
0.5142820-02 
0.4989600-02 
0.484517D-02 
0.470eeiQ-02 
0.457987D-02 
0.4457760-02 
0.4341950-02 
0 .4231980-02 
0.4127410-02 
0.4027  66D-C2 
0 .3932  98  D-02 
0.3642450-02 
0. 3755960-02 
0.3673300-02  ~ 
0.3594180-02 
0.3518380-02  ... 
0.3445700-02 
0.3375960-02 
0. 3308960-02 
0 .3244600-02 
0.3162630-02 
0,3123060-02 
0.3065640-02 
0. 3010290-02  . 

0.2956890-02 
0  .2905360-02 
0. 2855550-02 
0  .2807480-02 
0.276C98D-02 
0 ,2715980-02 
0. 2672430-02 
C. 2630240-02  .. 
0.258927D-02 
0.2549740-02 
0.2511310-02 
0. 2474020-02 
0.2437820-02 
0. 2402660-02 
0.236849D-02 
0.2335290-02 
0.2303000-02 


-C.5000C0D+00  ~ 
0.4440890-15 
0. 4085  21D-01 
0. 71 42  86D-01 
0.9141 69D-01 
0. 1046C0D+00 
0. 1 137320+00  - 
0.120424D+00 
0.125576D+0a 
0.1296910+00 
0. 1330640+00 
0. 1 35880D+00 
C. 1382630+00 
0. 1403C30+00 
0. 142067D+00 
0.1436050+00 
0. 144953D+00 
0.1461570+00 
0. 147223D+00 
0. 1481900+00 
0.1490600+00 
0.1499520+00 
0. 1505740+00 
0.1512360+00 
0. 1518460+00 
0. 1524090+00 
0. 1529300+00 
0.1534140+00 
C. 1538650+00 
0. 1542850+00 
0. 1546730+00 
0.1550450+00 
0. 1553910+00 
0.1557160+00 
0. 1560220+00 
0.1563100+00 
0. 1565  83D+00 
0.  1 56942D  +  00 
0.1570870+00 
C. 1573190+00 
0.1575410+00 
0.  1577510  +  00 
0.1579520+00 
0.  1581 44D+00 
0.1583270+00 
0.  1585C2D+00 
0.1506700+00 
0.  158831D+00 
0. 1 589E6D+00 
0.1591340+00 
C. 1592760+00 
C.  1594130+00 
0.  159545D+00 
0.  1596  72D+00 
0.  1597950  +  00 
0.1599130+00 
0.1600270+00 
C.  1601370+00 
0.1602430+00 
C. 1603460  +  00 
0. 1504450+00  . 
0.1605410+00 
0.  160634D  +  00 
0.1607240+00 
0.1609110+00 
C.  1609960  +  00 
0. 1609730+00 
0.1610570+00 
0. 1611350+00 
0.1612  10D+00 
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defined  in  section  2  of  this  paper. 
We  have: 


V 'a  i  4 


d  —  )  U  -  e  a, 


(47) 


;  where  in  (47)  expression  (17)  has  been  used,  and  as  in  section  4  the  conditioning 

X  >  1  is  deleted  from  expectations  involving  Y  .  Y  is  defined  in  section  2. 
q  —  3  SL 


A  simple  transformation  on  the  expectation  E 
equation: 


■fel- 


X  E{Y  }  +  E( 


XY 


XY 


leads  to  the  following 


(48) 


Xx 


Using,  as  in  section  4,  the  convexity  of  the  function  -  for  x  1,  and 

eAX-l 

recalling  that  Y ^  >  1,  we  obtain  the  following  inequality  from  (48) : 


XY 


.  -XY 
1  -  e  a 


<  X  E(Y  }  +  -4 

—  a  A 


(49) 


e  -1 


Substitution  of  (49)  in  (47)  gives: 


E<"/X  >!j  ^XE{V  +  ^ 


(50) 


Similarly,  using  the  concavity  of  the  logarithmic  function,  and  the  inequality 
(50) ,  we  obtain: 


!  "7 1 4  ^  e  { “A  i  4 1  e{y*} + 7Z ) 


(51) 


Now  recalling  that  is  the  same  with  Yd  in  section  2  of  this  paper,  and  using 
the  bound  (23)  in  section  4,  we  obtain: 
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At 

Finally,  using  expressions  (16)  and  (27)  in  section  4,  we  compute: 

•f‘7,  7 


(52) 


(53) 


Through  expressions  (46) ,  (50) ,  (51) ,  (52) ,  and  (53) ,  we  obtain  the  following 
upper  bound  cn  the  expectation  E-"m^},  where  the  conditioning  >_  1  is  now  deleted 
as  implicitly  necessary: 


(54) 


Due  to  difficulties  in  obtaining  lower  bounds  for  the  expected  value  of  the 
lower  bound  in  (46),  we  will  not  search  for  a  lower  bound  on  the  expectation  E{tn^}. 
Instead,  we  will  use  the  lower  bound  in  (29)  for  the  expected  value  e(d}.  As  we 
will  see,  our  numerical  results  justify  this  approach. 

The  e{d}  and  E(m^}  expected  values,  are  clearly  related  through  the  following 
equation: 

E{D}  =  .5  E{Ya>  +  E^i  (55) 

Substitution  of  the  bound  in  (54),  in  equation  (55),  provides  the  following 


upper  bound  on  the  expectation  E{d): 


Finally,  using  the  bounds  on  E{Y  },  as  given  by  expressions  (11)  and  (12)  in 

3 

section  3  of  this  paper,  we  obtain  from  (56)  an  upper  bound  on  the  expectation  e{d}, 
which  is  a  function  of  the  rate  X  only. 


In  table  3,  we  exhibit  our  computational  results,  where  the  upper  bound  is 
computed  through  expression  (56) ,  while  the  lower  bound  is  computed  through  expression 
(29)  and  is  the  same  as  in  table  1,  section  4.  Comparing  the  upper  bounds  in  tables 
1  and  3,  we  see  that  the  bounds  in  table  3  are  much  tighter.  In  addition,  we  observe 
that  for  high  X  values,  close  to  the  limit  .3465  for  stability  of  the  CCRA,  the 
lower  and  upper  bounds  in  table  3  are  very  close  to  each  other.  This  result  justifies 
the  fact  that  we  maintained  the  same  lower  bound,  as  in  section  4. 


0.0500 

0.100 

0.1500 

0.1696 

0.2000 

0.2500 

0.3000 

0.3333 

0.3400 

0.3450 


Lower  Bound 

Upper  Bound 

0.569 

0.615 

0.701 

0.814 

0.961 

1.193 

1.125 

1.431 

1.505 

1.979 

3.958 

5.266 

9.521 

11.889 

38.508 

42.157 

80.643 

85.119 

372.351 

380.126 

Upper  Bound 
From  Table  1 


15.096 

58.781 

122.043 

560.032 
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7.  Appendix 


Alternative  approach  to  the  computation  of 
Define: 

a  :  The  expected  delay  in  the  transmission  of  the  jth  successfully  trans- 

jN 

mitted  packet  within  some  CRI,  where  N  packets  are  involved  in  the 
collision  within  the  first  slot  of  the  CRI. 

Clearly,  a^  =  0. 

Now,  let  N  >  2,  and  let  the  first  trial  event  be  (i,  N-i) .  Then, 


jN 


1  +  a. .  ;  for  j  <  i  <  N 

Ji  -  ~ 


1  +  Li  +  5  0£i<j<N 


(A.l) 


Since  the  probability  with  which  the  event  (i,  N-i)  occurs  is  (!?)  — ^  ,  we 
obtain  from  (A.l): 


-  Z)  (V  a  +  Li  +  Vi.N-i’  +  2  (i>  (1  +  V 

i=j 


i-0 


i=0  i=0  1  3 


(A.  2) 


Applying  some  transformations  on  (A. 2),  we  finally  obtain  the  following  recursive 
expressions  for  the  a^k  parameters: 
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311  "  °*  a12  2*5,  a22  =  3,5 

N 

C2N-2)aiN  -  2N  +  1  +£  ( J)  «n  ;  N  >  2 

i=l 

j-1 

N-l 

*♦£(;) 

i=0  1=1 

i=j 

N-l  N-l 

-2“+Z)fi)Li+Z)<i)aii  : 

i=0  i=l 

N  >  2 

(A.3) 

2  <  j  <  N-l 
N  >  2 


Each  of  Che  N  packets  is  Che  ith  successfully  transmitted  packet  within  the 
CRI,  with  equal  probabilities  for  all  i  :  1  i  N.  Therefore, 


“n  = 


3  -VT 

lN 


(A. 4) 


By  substituting  expressions  (A. 3)  in  (A. 4)  and  applying  some  transformations. 


we  obtain  the  expressions  in  (37). 


